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Due to the internal structure of the nucleon, we should expect, in general, 
that the effective meson nucleon parameters may change in nuclear medium. 
We study such changes by using a chiral confining model of the nucleon. We 
use density-dependent masses for all mesons except the pion. Within a Dirac- 
Brueckner analysis, based on the relativistic covariant structure of the NN 
amplitude, we show that the effect of such a density dependence in the NN 
interaction on the saturation properties of nuclear matter, while not large, is 
quite significant. Due to the density dependence of the g^NN, as predicted by 
the chiral confining model, we find, in particular, a looping behavior of the 
binding energy at saturation as a function of the saturation density. A simple 
model is described, which exhibits looping and which is shown to be mainly 
caused by the presence of a peak in the density dependence of the medium 
modified aN coupling constant at low density. 

The effect of density dependence of the coupling constants and the meson 
masses tends to improve the results for E/A and density of nuclear matter 
at saturation. From the present study we see that the relationship between 
binding energy and saturation density may not be as universal as found in 
nonrelativistic studies and that more model dependence is exhibited once 
medium modifications of the basic nuclear interactions are considered. 



PACS numbers: 21.65.+f, 14.20.-c 



I. INTRODUCTION 

Since a nucleon is not a point object, but has structure, it must undergo changes when 
placed inside a nucleus. Among other properties, the meson-nucleon coupling constants 
may change. If this happens, it should affect the NN force. These effects have to be small. 
Otherwise, traditional nuclear physics, where one uses free-space two-body force would have 
failed badly. But even small changes in the NN force may have noticeable effect on the 
properties of nuclear matter. The purpose of this paper is to investigate possible changes 
of meson nucleon coupling constants due to the quark structure of the nucleon and their 
ultimate effects on the saturation properties, such as the density po and the binding energy 
—E/A of nuclear matter. 

The NN force in nuclear matter may become density dependent due to a wide variety of 
reasons. Any time one eliminates some degrees of freedom the resulting effective interaction 
becomes density dependent, the Brueckner G-matrix being the most widely known example. 
Another recent example is the work of Li et al on the effective interaction to be used in a 
mean field calculation which reproduces the results of a Brueckner-Hartree-Fock calculation. 

The density dependence studied here involves excitations of N* (/ = J = 1/2 resonances) 
degrees of freedom.^] Ultimately we are interested in the change of the NN force in the 
medium. In relativistic nuclear physics the most important forces are mediated by a scalar- 
isoscalar field a with rrio- ~ 600 MeV and a vector-isoscalar field uj with m^^ = 783 MeV. 
Thus the force range is short and the neighboring nucleons which can alter the internal 
structure of the interacting nucleons must also be close. Due to the Exclusion principle 

^The nuclear matter is an isoscalar and isotropic medium. It cannot change the spin or the 
isospin of a nucleon. It is also translationally invariant. But, as we will see later, it can still 
produce internal excitations in a nucleon. 
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and short-range correlations, the effective density of the polarizing nucleons, denoted as ^p, 
is less than half the normal nuclear density. It is quite possible that the A^* excitations 
produced by the neighboring nucleons may be treated perturbatively. However, an earlier 
study of this problem established that a fairly large number of resonances contribute. 
The mean field approach, which generates the eigen- combination of N and A^*s as the lowest 
state in the field due to the neighboring nucleons is a more expeditious way of calculating 
the effect. 

Our quantitative studies are based on the following strategy. We describe the structure 
of the nucleon with a model called the chiral confining model (COM) Specifically, we 

use the TOY Q version of this model. The role of nuclear matter is simulated with baths of 
external a and uj fields, the vacuum values being {(j).aac = ~F-k = —93 MeV and {uj)^ac = 0, 
respectively. The nucleon structure problem is solved in the presence of these bath fields. 
Then various properties of the nucleon, including meson nucleon coupling constants, are 
calculated for ranges of values of the two bath fields. 

The fact that the o and the uj fields are, in turn, produced by the nucleons through 
field- dependent coupling constants allows us to obtain density dependences of the o and 
the lo fields by solving appropriate nonlinear self-consistency equations. Once this has been 
done, we can use the known field dependences of various properties of the nucleon to obtain 
their density dependences. 

We indicate the density-dependent values of various properties of the nucleon with a star. 
Thus QaNN denotes the free space o NN coupling constant and 5'*7VAr the same quantity in 
nuclear matter. The latter is always dependent on p, the nuclear density. We may note that 
the properties of nuclear matter depend principally on glj^^, gl,^^, gl^j^, 5f*^^ and /*7VAr, 
where the last is the p Pauli coupling coefficient. 

A necessary input is the field-dependence of the exchanged meson masse^. We use the 

^The possible importance of density dependence of meson masses was first stressed by Brown and 
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simple model that the meson masses are linearly dependent on the a field, keeping the pion 
mass fixed. Q The coefficient of linear dependence of a and uj masses on the a field is chosen 
to give m* /m = 0.92 at normal nuclear density. 

The changes, as found in the CCM, clearly will have effect on the properties of nuclear 
matter. To study this we adopt the relativistic many-body approach [^|-|^ . Specifically, 
we carry out a relativistic Dirac-Brueckner calculation of the properties of nuclear matter 
using the one-boson-exchange model of Ref . . For the quasi-potential version of this model 



the full Dirac structure of the NN amplitude in free space has been constructed The 
resulting so-called 1A2 representation can be used to determine the saturation properties 
of nuclear matter P). Modifying the free space T-matrix to also include Pauli-blo eking 
and introducing our density-dependent meson-nucleon coupling constants, self-consistent 
relativistic Dirac-Brueckner calculations were performed in the manner of Ref. P] for a 
range of the effective density of the polarizing nucleons. 
Our main results are that 

(i) the effects of the density-dependent meson nucleon coupling constants, arising out of 
the quark structure of the nucleons, on the saturation density and —E/A at saturation 
are small but not negligible, 

(ii) and they do tend to improve the results. 

The next section contain brief introduction to the Chiral Confining Model (CCM) and its 
TOY version. Section 3 describes the calculation of nucleon properties as functions of bath 
a and u fields. These results are used in the next section to extract density dependences of 
nucleon properties. Section 5 describes the relativistic treatment of nuclear matter using the 



his collaborators ||5[ 

^Because of charge conjugation symmetry the masses of mesons, considered here, cannot depend 
linearly on uj . The dependence must be quadratic or of higher even power. 
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Dirac-Brueckner approach. The last section presents the main results, discussions of these 
results and concluding remarks. 



II. THE CHIRAL CONFINING MODEL 

An early attempt at extracting density dependence of nucleon properties using CCM is 
described in Ref. The present paper contains two significant improvements - inclusions 
of the instanton induced interaction and the pion cloud contributions to Qo-nn, 9pNN and 
fpNN coupling constants. The CCM has been described in detail in pj-Q. Here we give a 
brief review. 

The CCM is based on the notion of color dielectric function as introduced by Nielsen 



and Patkos [|rT|. By considering the average of all possible link operators, starting from 
X — e and ending on x with the paths completely contained in a four-dimensional hypercube 
of side L, they introduced a color singlet, Lorentz scalar quantity K{x) and a color octet. 



Lorentz vector, coarse grained gluon field 5^: 



K(x) = lim — Tr 



1 d 
-yXaB'' = hmi— 



Upon integrating out the QCD gluon fields in favor of these new collective variables one 
obtains the Nielsen-Patkos lagrangian in the form of a derivative expansion: 

C^P = ^{x)[zK{x)^ 5 -K{x)m^ - gl^{xmx) - ^G^^G'^'^'' + .... (1) 

The gauge field is B^/ K and not B^. The coarse grained field tensor is 

Tja na Tjb r>c 

G% = d,^-d.^ + r'^^^. (2) 

From the gluonic term one identifies e = K"^ as the color dielectric function. Nielsen and 
Patkos conjectured that 

{K),ac = 0. (3) 
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This conjecture, crucial for our model, has been justified from the lattice gauge point of 



view by Lee et al |12 



A quark has ever-present interaction with the quark condensate of the vacuum. If 
{K)j,ac = 0, the interaction will appear to be infinitely strong compared to quark kinetic 
energy and a quark cannot exist in that region. It can only reside in the region where 
{K)yac 7^ 0. A color singlet quark system can polarize the vacuum and change the value of 
K away from zero, thus dynamically generating the bag where the quarks can stay. 

The Nielsen-Patkos lagrangian recognizes the existence of gluon condensate through the 
vanishing of {K)vac- However, the quark condensate is not manifest. Without it one cannot 
develop an effective lagrangian which contains the physics of the interaction of a quark with 
the quark condensate. We deal with this problem by conjecturing that one can integrate out 
the coarse grained gluon fields in favor of meson fields as new collective variables P, p!3| , p!4 



It is also necessary to introduce, following Nielsen and Patkos |jTT|, a new color singlet, 
Lorentz scalar field, Xi proportional to K by the equation 

K{x) = g^x{x). (4) 

Being related to K, a purely gluonic object, the x field is a member of the glueball family. 
Hence, it is a chiral singlet, a fact evident from the first term of the Nielsen-Patkos lagrangian 
given by Eq. (|l|). Large N^oior analysis shows that it is a hybrid field @,|15|. The result 
imposes the requirement that 

{X)vac = 0. (5) 

Retaining minimum powers of fields and their derivatives the basic lagrangian of the 
CCM has the form: 

1 ^ 

CccM = K{x)^{x)[i- ^ -mg]^{x) 

+ C„ + ^d^xd^X-U{x). (6) 



The quantities (?7r, Qui and Qp are the quark- meson couphng constants. The quantity niq is 
the current quark mass. Its value is set at 7.5 MeV. It contributes 17 MeV to the vrN sigma 
term Because of its negligible role in the present work we will neither refer to this term 
nor count it as a parameter in our subsequent discussions. However, it is included in the 
actual numerical work. 

We use the chiral invariant lagrangian of Lee and Nieh [T^ for Cmeson- The lagrangian 



ensures that {cr)yac = —Ft^ and that the mesons have their respective physical masses when 
calculated at the classical ( or tree) level. The fields vr and a form a (1/2 x 1/2) representation 
of SUl{2) X SUr{2), while the fields p± Ai form (1,0) and (0, 1) representations. 

To complete the definition of Cqcm one must specify the x potential. Since nothing 
substantial is known about it, we try two simple forms: 

1 



PUREMASS: f/(x) = -m^x', 

1 



QUARTIC : U{x) = ^^xX'(l - x/Xof- (7) 

The QUARTIC potential has two minima, the one at x = describes the vacuum while 
the other at x = Xo is an 'excited' state, which for simplicity we keep degenerate with the 
vacuum. The value of xo, the location of the second minimum is set at 40 MeV. The hybrid 
mass, is set at 1400 MeV. As we will discuss later, the results of the calculations depend 
largely on one particular combination of these parameters and rather weakly on individual 
ones. 

Any mean field calculation containing isovector fields requires using states which do 
not have good isospin symmetry. To accomplish this we use hedgehog spinors and fields, 



introduced first for tt and cr fields by Chodos and Thorn and extended to vector meson 
fields by Broniowski et al [|18[. We closely follow the mean field analysis of Ren et al ||19|| . 

It should be emphasized that the appearence of K in the denominator of the quark- 
meson interaction term is not ad hoc. It is obtained by matching the K dependence of the 
four-quark interactions which arise on the one hand from C^phy integrating out the 5^ 
fields and on the other hand from Cccm by integrating out the meson fields P, |20|j2T[| . 



The presence of the factor K with ipi}^ ^ ip changes the field canonically conjugate to 
%l) from the usual iip'^ to iKip'^ . As a result the quark term in all Noether's currents carries 
the factor K. The transformation \/KiI){x) ip{x) makes the pair ijj and iip"^ canonical 
and removes the factor K from all Noether's currents. Two changes occur in Cqcm- The 
free quark term no longer has the factor K, while the quark meson interaction acquires the 
factor in place of K in the denominator. 

Mean field treatment of Cccm with a variety of reasonable sets of parameters re- 
veals that even where the quark density is large the meson fields differ only slightly 



from their respective vacuum values. This suggests strongly that we introduce a simplified 
version of CCM. The simplified version, which we call the TOY model, consists of fixing all 
meson fields at their vacuum values: 

{T^)vac = {^)vac = {p) vac = (^l)mc = 0, (8) 

leaving only the quark fields and the x fi^ld as dynamical variables. The Toy model la- 
grangian, in terms of the canonical quark field, is given below. 

^Toy = ^{x)[^l 5 -m, - ^^W^) + \d,xd'x - U{x). (9) 

Not counting m^, the TOY model has only two parameters, Qt^F-^/q'^ and in the PURE 
MASS version and one additional parameter, xo, in the QUARTIC version. Formation of 
the bag is even more transparent in the TOY model as the constituent quark mass term, 
JjjjrF}^, becomes infinite when v 0. 

The unusual form of the quark- meson interaction with the factor of = {QxXY i'^ ^^e 
denominator is of some importance in the present investigation. In its absence a constant 
external uj field will merely shift the energy of the quark without exciting it to higher 
states. Thus the nucleon will not be polarized. A constant a field will excite the quark. 
Orthonormality of different eigenstates of the Dirac hamiltonian ensures that / dP'ru'^^u = 0, 
but not the vanishing of / d?ru'u. The latter gives the effect of a constant external a field. 
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The presence of the factor {QxX)'^ enables a constant external uj field to excite the quark and 
enhances the ability of a constant external a field to do the same. 

Because of its importance it is necessary to provide some evidence for the presence 
of iOxX)'^ ill the denominator of the quark-meson interaction term. Fortunately, there is 
a verifiable consequence of this feature. Gauge invariance ensures that there is no such 
modification of the photon-quark interaction. Thus contrary to the conjectures of universal 
coupling or current-field identity |^ the isovector couplings of the photon and the 
p meson differ by this factor. Because of this difference, CCM predicts that the p- 
nucleon Pauli coupling constant, fpNN should differ from the isovector anomalous magnetic 
moment, ni=i. The predicted [|T9| value is fpNN/'^i=i = 1-4, instead of 1 predicted by 
the two conjectures. Based on dispersion theoretic analysis of vrA^ scattering data Hohler 



and Pietarinen ||2^ had estimated the ratio to be ~ 1.75. Our value is within the possible 
uncertainties of this result p5[|.^ 

The role of the pion field on baryon properties is never negligible. In the TOY model one 
takes account of it perturbatively P,pB[ wherever appropriate. In the context of the present 
study such corrections are important in the calculations of the coupling constants of mesons 
of even G-parity, gaNN, QpNN and fpNNf^M. 

The Toy model with perturbative pionic correction gives poor results for the masses of 



and A and for the nucleon charge radii 0j2^. McGovern et al [|23 showed that there 



is approximate, but very good scaling behavior in the TOY model. If one chooses the pa- 
rameters to get the nucleon mass right the values of most other quantities are nearly fixed. 
Thus the bad results for mass and size cannot be improved simultaneously by varying pa- 



^ In the MIT bag model also the photon and the p meson couple to the quark in different manner. 
Brown and Rho had exploited this feature to explain the vahie of fpNN/f^i=i- Dialing the chiral 
angle they fitted the value 1.75. Our result is approximately independent of parameters because 
of the scaling property mentioned later. 
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rameters. One needs a new interaction which is short-ranged and more attractive between 
singlet-singlet quark pairs than between triplet-triplet pairs. The magnetic one-gluon ex- 
change interaction ||2^ fits the requirement exactly. However, the physics of exchange of a 
color octet tower of gluons is already incorporated in meson exchanges. We cannot include 
it again. 



We solve the problem by including the instantanton induced 't Hooft interaction ||29| , |30 
In the two fiavor case of u and d quarks, the 't Hooft interaction generates interactions 
between quark pairs in fiavor antisymmetric state only and the interaction is attractive and 
zero-ranged. So it fits our requirement without duplicating any mechanism already included. 
Using Shuryak's ||31| description of the vacuum as an instanton liquid the strength of the 't 



Hooft interaction can be estimated in QCD, but not in CCM 0. So we have to treat it as 
a free parameter. We fix it by fitting the — A mass splitting with the 't Hooft interaction 
together with one-pion exchange contribution 0. 

The method of calculation with the TOY model in the absence of external a and lj 
has been described in detail in Ref. and will not be repeated here. The changes due to 
the presence of external fields are straightforward 0. Going back to Cccm, described by 
Eq. (j^), we interpret a and Uq as the bath fields. The vacuum expectation values of all other 
fields, TT, Ui [i = 1,2,3), and Ai^ are zero because nuclear matter has good isospin and 
good parity. The quark field and the x fi^ld continue to be the only dynamical variables 
in the extended TOY model lagrangian. The new energy functional of a nucleon in a a -a; 
bath is given below. 



E = N^ I d\u\r) 
■1 



-la ■ V - 



u{f) 



{9xX{r)y 

+ J d^r 

- 2CsJ d^rZ{r){G\r) + F\r)). (10) 



where the valence spinor, M(r) is 
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c, 



(11) 



G{r) 
^ id ■ fF{r) ^ 

and C is a 2-component Pauli spinor. 

The last line in Eq. ( p!OD represents the contribution of the instanton induced 't Hooft 
interaction with a regulating function, Z{r), which cuts off the interaction at a distance 
scale of 0.25 fm. The details may be found in Refs. PJ26[1. 

Note that the external fields are constant in space and time. The vacuum value of the a 
field is —Ft,. Nuclear matter generates a field with positive sign thus reducing the magnitude 
of the net a . It is customary to represent the effect as a modification of the pion decay 
constant: 



^ F-ix ~\~ (^nmatter . 



(12) 



The uj in Eq. (|ToD is the time component of vector uj field. The static nuclear matter 
distribution can generate only the time component. However, a nucleon moving with velocity 



V in nuclear matter does see a space component Q = —vu/y/l — v^. We ignore the role of 
this term. 



III. FIELD DEPENDENCE OF NUCLEON PROPERTIES 

We pick values of F* in the range 63 MeV to 93 MeV in steps of 6 MeV and values of lo 
in the range to 40 MeV in steps of 8 MeV. For each pair of F* and u we find the spinor and 
the X fisld which will make the energy functional stationary. The condition of stationarity 
yields coupled nonlinear equations. Ref. describes the method of solving these equations. 
Once the solutions are obtained we can calculate the desired properties of the nucleon in 
the a -UJ field bath. 

We introduce the dimensionless variables x and y to describe, in units of m^, the bath 
fields a and uj : 
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(y - {(y)vac ^F^-F*^ xm^. (13) 

LO = ym^, (14) 

We fit the results of every pliysical quantity Q of interest witli tlie quadratic form: 

Q{x, y) = a[l + hx + cy + dx^ + fxy + gy"^] = aF^(x, y). (15) 



TABLE I. Values of quantities a, b, c, d, f and g for the PUREMASS case. The quantities 
under the column a are the free space values. 



Quantity 


a 


b 


c 


d 


f 


g 


{g^NN/2My 


0.83 


1.72 


-0.78 


3.51 


-5.71 


0.13 


qIjNN 


11.02 


1.10 


-1.60 


2.44 


-2.90 


3.22 


(q)* 
9aNN 


2.34 


3.14 


-5.25 


3.34 


-21.77 


6.99 



TABLE II. Values of quantities a, b, c, d, f and g for the QUARTIC case. The quantities under 
the column a are the free space values. 



Quantity 


a 


b 


c 


d 


f 


g 


{g^NN/2My 


0.92 


1.62 


-0.49 


3.55 


-3.36 


0.045 




13.72 


0.74 


-0.93 


1.23 


-1.18 


1.45 


9aNN 


4.92 


1.67 


-2.13 


-2.39 


-6.65 


1.72 
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The results of our fit are shown in Tables | and for the PUREMASS and QUARTIC 
case respectively. The quantity a is the value of Q in vacuum and F^{x, y) is the quadratic 
polynomial in square brackets. The three coupling constants of special interest are g^^Jf^, 
9l)NN {qt^nn I'^M)* . The a-nucleon coupling constant carries a superscript (g) to indicate 
that only the direct coupling of the cr field to valence quarks is included in glj'}j*N- The full 
coupling constant also includes coupling to the pion field produced by the valence quarks. 
The u and the vr fields have odd G-parity. Hence their source currents couple to 3 or higher 
odd powers of the vr field. Since the vr field is weak inside the nucleon, the higher power 
contributions are not important and have been omitted in the present calculation. The 
superscript (g) is redundant for g^j^j^ and {gT^^j^f /2M)* and have been omitted. 

The parts of the meson source current which show the coupling to quarks may be read 
off Eq. (H) describing Cqcm- We list them below: 

ii'>=9.^. (16) 

(17) 

.S^..^. (18) 

The meson-nucleon couphng constants are defined in the usual way. The expressions listed 
below are for p — 0. 

{proton((p)) tl j^p'f^ 



{prot(m{{p)) I d^f [rxjj^h 



1 Niip))) 




(20) 


1 Niip))) 


= gujNN, 


(21) 


proton{{p)) 1) 


1 (g) 


(22) 


proton{{p)) 1) 




(23) 


proton{{p)) 1) 


= {p -p)zg^NN/2M. 


(24) 



2 

{protonCip)') Tl jS(0) 
The quantity HpNN appearing in Eq. (^) is the p- magnetic moment of the nucleon: 

^J^pNN = (gpNN + fpNN)/'^M. (25) 
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We list below the results from CCM for the contributions to the various coupling con- 



stants coming from direct coupling of meson fields to the quarks. 

gitN = 9^1 d'r[G'ir) - F\r)]/K\r\ (26) 
1 n r 

ig^NN/2My = g^- J d'rrG{r)F{r)/K\r), (27) 

glNN = 3^7. / d'r[G\r) + F\r)]/K\r), (28) 

9% = 9,j d\[G\r) + F\r)]/K\r\ (29) 

= 9PY I dhrG{r)F{r)/K\). (30) 



To obtain density-dependent, or equivalently field-dependent, coupling constants one must 
use K and spinor components G and F obtained in the presence of the bath fields. The 
stars should be removed for free space values of the coupling constants obtained with G 
and F appropriate for free space. 

The nucleon mass term which appears as a divisor for Qt^nn and for the p magnetic 
coupling stands for the experimental value of the mass, i.e., 939 MeV. It is used to define 
dimensionless coupling constants from the integrals of dimensions of length which appear in 
Eqs. ( P7D and (PU]). Our studies measures the density dependence of the quantities repre- 
sented by the integrals. 

A model of hadrons cannot anticipate what effective hadronic lagrangian its results will 
be used for. Hence, to remove possible misunderstanding we have stated in Eq. (plf ) the 
precise definition of the quantity gj^NN/^M as it appears in our work. The expression on 
the left hand side involves the pion source current, usually well-defined in most quark based 
model of the nucleon. It is also well-defined in any effective hadronic lagrangian. If the 
effective lagrangian uses pseudoscalar vrN coupling the quantity g-^NN is the dimensionless 
coupling constant of the theory. If, on the other hand, it uses derivative coupling then 
g-TrNNf^M is the coupling constant of the theory, having the dimensions of length. Note, 
again, that M in the denominator of (7^7VAf/2M for the derivative coupling theory is a fixed 
number by definition, usually 939 MeV. 
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The mesons with even G-parity, a and p, also couple to the pion cloud, the contributions 
from which depend essentially quadratically on ((^jtatat/SM)*. Hence we write: 

qInn = g^L + U9^NN/2Mr\ (31) 
qInn = 9% + U9.NN/2Mr\ (32) 

^^;NN = f^%N + Cfi9nNN/2My'. (33) 

The parameters, (g and C/, represent the strengths of the pion cloud contributions. The 
quantities and (g have the dimensions of mass squared , while (f has the dimensions 
of mass. We choose reasonable values of g^NN, 9pNN and fpNN/^M in free space and use 
the results of the TOY model in free space for S'^tvat' 9'pNN^ l^^pNN ^'^d gT,NN/2M to fix the 
coefficients (g and (f. These values and the relevant parameters are listed in either 



Table |II| or in the table caption. 

We will see later in section that the free space meson nucleon coupling constants used 
here are not exactly the same as in the one boson exchange potential used in the Dirac- 
Brueckner calculations. However, the differences are small and are not expected to affect 
the qualitative purpose of this paper. 

As stated in the introduction a necessary input is the field-dependence of the exchanged 
meson masses. We use the simple model that the meson masses are linearly dependent on x 
1^, keeping the pion mass fixed. We assign the same x-dependence to both a and uj masses. 

m* = m(l - bmesonx) = m(j){x). (34) 

The function 0(a;) has been introduced so that the rest of the discussion is not specific to 
the simple linear dependence used in the numerical work. The coefficient bmeson is chosen to 
give m*/m = 0.92 at normal nuclear density. Its values are listed in Table |T|. 
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TABLE III. List of parameters and coefficients. Both the PUREMASS and QUARTIC cases 
use nia = 650 MeV, niu, = 783 MeV, nip = 770 MeV, m-„ = 140 MeV, g^NN = 8.0, gpNN = 5.0, 
l-i'pNN = 2.45 and m*/m = 0.92 at normal nuclear density. The numbers in the table are all in 
pion mass units. The quantity ^ttAtat/SM has the dimensions of length and is expressed in units 
of m~^. The quantities g^jNN, qI^nn^ 9pNN ^"^^ bmeson are dimensionless. The quantity Q has the 
dimensions of mass and is expressed in units of m^. The quantities and Q have dimensions of 
mass squared and are expressed in units of m^. 



QuMiiiitA- 


PITREMASS 


QITAETIC 


g^NN/2M 


0.8349 


0.9241 


9luNN 


11.02 


13.72 


M 

9aNN 


2.389 


4.919 




8.049 


3.608 


9^pNN 


2.448 


3.048 




3.661 


2.286 


t^pNN 


0.8349 


0.9241 


Cf 


2.317 


1.787 


bmeson 


0.4769 


0.4622 



16 



IV. DENSITY DEPENDENCE 



The first step in converting the information about field dependence of the couphng con- 
stants to density dependence is to obtain the density dependence of the fields, a and u; 
, or equivalently, x and themselves. This is done through the following two basic self- 
consistency equations: 



X 



9aNN 



m 



*2 



y = —^^P' 



(35) 
(36) 



where the various quantities have been defined by Eqs. (|T5|) and (|3^). 

The factor ^ represents the reduction of the effective density of nucleon available to 
polarize a given interacting pair. The reduction is the combined effect of the Exclusion 
principle and short-range correlations among nucleons. An examination of the illustrative 
Goldstone graphs of Fig. |l| may be useful to understand the origin of this factor. 




FIG. 1. The Fig. (a) represents interaction between two Fermi sea nucleons interacting via a 
G-matrix. In Fig. (b) a third nucleon excites one of the interacting nucleons into a resonance 
state, leaving a hole in the Fermi sea. The resonance propagator is shown as an upward line with a 
star. Fig. (c) is the result of exchange between the third nucleon and the left interacting nucleon. 
Recall the convention that a line starting and ending at exactly the same time represents a hole 
line. Fig. (d) represents the exchange between the third nucleon and the right interacting nucleon. 
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Besides the diagrams shown in the figure there are others obtained by interchanging the 
roles of the right and the left interacting nucleons. There are also diagrams where a fourth 
nucleon excites the interacting nucleons to a resonance state, generating two resonance-two 
hole states. The two resonance states interact and fall back into the hole states. All these are 
examples of single interaction of the neighboring nucleons with the interacting pair. These 
may be repeated and summed to generate field effects in a space containing the nucleon and 
J = J = 1/2 resonance states as degrees of freedom. If resonances were not included the 
procedure will generate contributions to the usual average field field of a many body theory 
where only the nucleon degree of freedom is counted. The mean field of CCM describes 
approximately the multidimensional average field described here. 

Inclusion of all exchange diagrams enforces the Exclusion principle, which ensures that 
only two other nucleons can come close to the two interacting nucleons. The latter are 
mostly within a short distance. This consideration alone makes ^ = 1/2. Repulsive short 
range correlations make ^ even smaller. We expect that < < 0.5. A full self-consistent 
many body treatment should determine this factor. But for the present preliminary study 
we take it to be an unknown parameter confined to the range < ^ < 0.5. 

In Eq. ( P^D one should use the scalar density. This can be done only through a full many 
body self-consistent calculation. We simplify our work by using the vector density p. This 
simplifying strategy allows us to factor the problem of extracting the density dependence of 
the coupling constants from the rest of the many body problem. We estimate the error to 
be ~ (p^)/4M^ ~ 0.05. The strategy is probably justifiable in a preliminary study of the 
main question studied in this paper. 

The Eqs. (|35D and (^) may be rewritten as: 



X = 



{S('^)F'^(x,i/) + S-(F-(x,2/))2} 




(37) 



where the combinations S'^'^^ and S'^ are 



9(tNN 



'-fpo and = C.((7.iV7v/2M) 




(38) 



18 



The quantity po is the normal nuclear density. The quantities T,'^ and S'^ would be the 
a fields produced by nuclear matter due to direct a -quark coupling and due to a -pion 
cloud coupling, respectively, if the coupling constants did not change with density. The x, 
y dependent quantities were defined by Eqs. (|1^) and (Q). 

For future convenience we introduce the dimensionless effective density p: 

P = e-. (39) 

Po 

In a similar manner we write 
where 

^ = — 5-P0) (41) 
mi 

would be the u field produced by nuclear matter if the coupling constant did not change 
with density. 

We solve the two coupled nonlinear Eqs. ( |T7D and (^D|) for a range of p giving us the fields 
X and y as functions of p. Using the results and Eqs. (|T3|) we obtain numerical values of the 
five coupling constants, g*NN^ 91nn^ {g^rNN I'^M)* , g*p^N, p*^^ and m*/m as functions of p. 

We remind the reader that we do not use the absolute coupling constants from the CCM 
calculation in the Dirac-Brueckner treatment of nuclear matter. We use only the results for 
the density dependences of the ratios such as qInn / Q^nn- 

These ratios are plotted as functions of p = ^p/po in Fig. 0. To make it convenient to 
use these results in a relativisitic Dirac-Brueckner-Hartree-Fock calculation we represent the 
density dependences in the form of a [4, 5] rational function: 



9 1 + E„=i,5/3^P"' ^ ^ 



m l + En=l,5/?rP"' 
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The coefficients «£ and /3„ for the five couphng constants and the meson mass ratio are 
obtained by least square fitting and are hsted in Tables ^ and [V|. 




1.20 
„ 1.15 

a 

"S 1.05 
g' 

i. 100 
g 

" 0.95 
0.90 



(c) 

QUARTIC 










tj 



.0 .1 .2 .3 



1.20 

» 1.15 
rl 

3 1.10 

o 

c 1.05 

"H. 

S 

u 1.00 
0.95 



(b) 








PURE MASS 








t 

/» . - 
















.1 


.2 .3 
fP//>„ 


.4 .5 


(d) 








QUARTIC 








f 

p 











.2 .3 
tp/p„ 



FIG. 2. Plots of g*^^]^ / QcrNN labelled with a , {gT^NN f^^)* / (dnNN /'2M) labelled with vr, 
m*/-"^, gljNNlduNN labelled with w, g*pNN/9pNN labelled with gp and HpNN/if^pNN) labelled with 
fp as functions of (,pn/po- Results for both the PUREMASS and the QUARTIC forms for C/(x) 



are presented. All relevant parameters are listed in Table III in section III 
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TABLE IV. Coefficients of rational functions defined in Eqs. (|4^) and obtained witli PURE 
MASS U{x). Columns 2 through 6 are for the ratios of coupling constants, while the last column 
is for m* /m. Rows labelled 1 through 4 are the numerator polynomial and the remaining rows are 
for the denominator polynomial. All parameters relevant to the calculations are listed in Table III 



in section III 





9crNN 
gtjNN 


guiNN 


(gnNNfiM)* 

g„NN/2M 


9pNN 
gpNN 


^^pNN 

MpJVJV 


m* /m 


1 


0.1553 


-0.0779 


0.1817 


0.1491 


0.3015 


-0.0877 


2 


-0.4452 


0.2303 


-0.4327 


-0.2797 


-0.7174 


0.1535 


3 


0.4843 


-0.2198 


0.4802 


0.2101 


0.7951 


-0.1481 


4 


-0.2853 


0.0484 


-0.2026 


-0.0700 


-0.3352 


0.0295 


5 


-2.3792 


-1.3216 


-2.6652 


-2.6050 


-2.7329 


-2.1338 


6 


2.8749 


0.3065 


3.6722 


3.1574 


3.8297 


2.6547 


7 


-1.7181 


1.0527 


-2.6344 


-1.8743 


-2.8004 


-1.6752 


8 


0.6688 


-0.6230 


1.2151 


0.5764 


1.2733 


0.9508 


9 


-0.0517 


0.1177 


-0.1483 


-0.0486 


-0.1415 


-0.1372 
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TABLE V. Coefficients of rational functions defined in Eqs. (H,® obtained with QUARTIC 
U{x)- Columns 2 through 6 are for the ratios of coupling constants, while the last column is for 
m* /m. Rows labelled 1 through 4 are the numerator polynomial and the remaining rows are for 



the denominator polynomial. All parameters relevant to the calculations are listed in Table III 



m 



section III 





9(tNN 

Qa-NN 


9ujNN 


g^NN/'iM 


9pNN 
QpNN 


t^pNN 
MpiVJV 


m* /m 


1 


0.0522 


-0.0653 


0.1928 


0.1117 


0.3129 


-0.0852 


2 


-0.2078 


0.1460 


-0.3308 


-0.2173 


-0.5382 


0.1093 


3 


0.2027 


-0.1240 


0.2994 


0.1304 


0.4855 


-0.0842 


4 


-0.1174 


0.0297 


-0.1087 


-0.0230 


-0.1756 


0.0140 


5 


-1.9719 


-1.0893 


-1.9685 


-2.7149 


-2.0459 


-1.5207 


6 


2.0516 


0.2868 


2.2234 


3.1836 


2.3486 


1.4465 


7 


-1.0666 


0.5555 


-1.2896 


-1.9596 


-1.3931 


-0.6018 


8 


0.3575 


-0.2999 


0.5296 


0.5796 


0.5573 


0.2730 


9 


-0.0292 


0.0451 


-0.0640 


-0.0663 


-0.0601 


-0.0200 
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V. DIRAC-BRUECKNER ANALYSIS 



Modifications of tlie meson masses and meson-nucleon coupling constants by the presence 
of the nuclear medium as discussed in the previous sections will, in general, lead to changes 
in the saturation properties of nuclear matter. This problem can be studied within the 
relativistic Dirac-Brueckner approach. To carry out such a study we follow the approach 
described in Ref. As dynamical input we use the relativistic one-boson-exchange model 
of Ref. P|. The NN T-matrix satisfies a Bethe-Salpeter equation, which is formally given by 



T = V + VS2T, 



(44) 



where V is the NN interaction and 5*2 is the free 2-nucleon Green function. In this study we 
use in particular the quasi-potential version of it. The interaction V is assumed to be given 
by the exchange of vr, a, p, uj, 6 and rj mesons. In particular, in the model a derivative pion 
and 7] meson coupling has been assumed. For our study we choose the coupling parameters 
of interaction A, which gives a reasonable fit to the NN scattering phase shifts To 
regulate the behaviour at high momenta, a monopole form factor A^/(A;^ — A^) has been 
introduced at each meson-nucleon vertex. A cutoff mass of A = 1150 MeV has been taken. 
In Table |V1| the meson parameters of the model A are listed. 

TABLE VL Meson-nucleon coupling constants g and meson masses m of the 
one-boson-exchange model A 





TT 


a 


UJ 




P'/P" 


6 


7] 




14.2 


7.6 


11.0 


0.43 


6.8 


0.75 


3.09 


m (MeV) 


139 


570 


783 


763 


763 


960 


548 
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According to the CCM model the free space meson parameters will be modified in the 
presence of the nuclear medium. Within the one-boson-exchange model the medium modifi- 
cations found in the previous section can be simply implemented by replacing for each meson 
(f) the coupling constant q^iNN and the mass m in the one-boson-exchange interaction V by 
the corresponding g^j^j^ and m* given by Eqs. and (^Sp. Moreover, the Pauli-blocking 



due to the medium has to be accounted for in Eq. (^4|), leading to the Bethe-Brueckner- 
Goldstone equations for the G-matrix. This is done by replacing the Green function in 
Eq. (H 

^2 ^ Sip,)Sip2)Q (45) 

where S{pn) are the medium-modified nucleon propagators and Q the Pauli-blocking opera- 
tor, which projects out in the intermediate states nucleon momenta inside the Fermi sphere. 
Introducing the baryonic current B = pu, with p being the density and u the unit vector 
u = (1,0) in the nuclear-matter frame, due to Lorentz covariance we may write the nucleon 
self-energy contribution in terms of the 3 invariants 

S(p) = - S°7.u - S''7.px, (46) 

where p± = p — {p.u)u. The medium-modified nucleon propagator can be approximated 
near the dressed nucleon pole as 

S{p) = [E*^o-Pl-M*]-' (47) 

with 

M* = (M + S')/(l + S^) (48) 
E* = (po + S°)/(l + S^) 

For the Pauli-blocking operator Q an angular averaged approximation has been made. The 
resulting relativisitic G-matrix equations 

G = V + VS{q,)S{q2)QG (49) 
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with Qn the momenta of the nucleons in the intermediate states, are solved in the A^A^ CM 
frame after partial wave decomposition. This is done within a quasi-potential equation de- 
scription using the helicity basis of positive and negative energy spinor states corresponding 
to mass M*. 

Relativistic nuclear matter calculations require the knowledge of the transformation 
properties of the G-matrix under Lorentz transformations. After having determined the 
relativistic G-matrix in the 2-particle CM system at a fixed matter density, we use the IA2 
representation [§] to obtain it in the nuclear matter frame. This representation gives the 
complete covariant form of the G-matrix in an unambiguous way. To reconstruct it all matrix 
elements of the amplitude in the full Dirac space are needed. From this IA2 representation, 
neglecting the vacuum fluctuation terms, the self-energy can be readily determined. Taking 
nucleon 2 to be one of the valence nucleons in the Fermi sphere with Fermi momentum pf, 
we have 



where nucleon 2 is in a positive energy state with momentum p2 and Tr2 designates the 



summation over the spin index of this nucleon and E2 = ^p\ + M*^. Hence the self-energy 

contribution is essentially obtained by taking the diagonal matrix-elements of the G-matrix 
with one of the nucleons belonging to the filled Fermi sphere while the other one carries a 
momentum p. 

Since the G-matrix depends implicitly on the S"'s through Eq. (^) the calculations have 
to be carried out in a self-consistent way. In doing so, the self-energy can be determined for 
a given matter density p. In particular we have assumed that S does not vary significantly 
within the Fermi sphere, so that its value can be taken at the Fermi momentum pf. Some 
studies 1^ of the momentum dependence have been made, indicating that the variations in 
S are of the order of 10 %. Hence this is a reasonable approximation. 

The partial wave decomposed integral equations for the G-matrix have been solved using 
Gaussian quadratures. Typically 24 Gaussian points are sufficient to get an accuracy of a 




(50) 
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few percent. Moreover, 6 Gaussian points for each integration variable in Eq. (^0]) have 
been used to determine the self-energies. The self-consistent solutions can be found in an 
iterative way. Adopting a starting value for the S's in Eq. (0), Eq. (^) is solved to yield 
through Eq. ( |50D new values for the S's. Varying subsequently the S's we then can determine 
the solutions of Eqs. (^^^D|) such that the S's are the same. Once the selfconsistent self- 
energy solutions have been found the binding energy of the ground state can readily be 
calculated from the energy-momentum tensor as a function of density. For more details of 
Dirac-Brueckner calculations we refer to Refs ^ and pf. 

For the various relativistic one boson interactions we find, in general, that the system 
exhibits saturation. As found in Ref an exception to this is when one assumes that only 
the meson-masses drop as m*/m = M*/M [Q. It is interesting to note that in the case we 
also allow for medium-modification of the meson coupling constants as is predicted in our 
study of the chiral confining model, saturation does occur, that is we find that the binding 
energy of the ground state as a function of the matter density p has a minimum. In the next 
section we discuss our results for the saturation properties of nuclear matter as predicted 
for the medium modifications of the coupling constants as found in our COM model. 



VI. RESULTS AND DISCUSSIONS 
A. Main Results 

As we have discussed earlier, the coupling constants are dependent not directly on the 
density p of nuclear matter, but on the effective density, ^p, seen by an interacting pair in 
its immediate neghborhood. The Pauli Exclusion principle ensures that ^ < 1/2. There is 
certainly further reduction of ^ due to strong repulsion at short distances. Thus the range 
of ^ of interest is < ^ < 0.5. 

We pick a value of ^ and proceed to calculate E/A and p at saturation using the Dirac- 
Brueckner approach described in Section 0. In the calculations reported here we have 
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used the one boson exchange interaction A specified in Ref. P]. During each loop of self- 
consistency iteration we use coupling constants and meson masses appropriate for p^r = 
C,p/po, defined earlier by Eq. (|39D, p being the nuclear density at the particular stage of 
the iteration. Specifically, we use Eqs. (^) and ( ^31) with the coefficients of the rational 
functions listed in Tables |rV| and |V|. Finally, self-consistency is achieved and we have a pair 
of numbers, E/A and density p at saturation for the given value of ^. Figures |^ and ^jshow 
these results for the PUREMASS and the QUARTIC cases, respectively. 

1 1 1 1 1 1 — 

PUREMASS J 

.0 - 

" ^ .1 



1 

■8 ^ - 



0.0 0.1 0.2 

P (fm ) 

FIG. 3. Results of E/A and p at saturation for the PUREMASS case for a series of values of 
< ^ < 1.0, shown with solid circles. The square box represents the empirical nuclear matter 
density and binding energy. The values of appear next to the solid circles. The line is for guiding 
the eye only. The details of calculations are given in section 0. Only the range < < 0.5 is of 
physical interest. 
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0.0 0.1 0.2 

P (fm ) 

FIG. 4. Results of EjA and /o at saturation for the QUARTIC case for a series of values of 
< ^ < 1.0, shown with solid circles. The square box represents the empirical nuclear matter 
density and binding energy. The values of ^ appear next to the solid circles. The dashed line is for 
guiding the eye only. The details of calculations are given in section Only the range < ^ < 0.5 
is of physical interest. 

The value corresponding to ^ = represent the Amorim and Tjon [§] result. From the 
figures we see that dependence of the nuclear saturation properties on ^ is just of the right 
size to be of interest in the present study. The E'/A vs. p curve for the PUREMASS case 
passes through the square representing experimental data for ^ slightly greater than 0.1. 
For the QUARTIC case the EjA'vs,. p curve comes close to the square, but does not enter 
it. The closest approach occurs for ^ a little larger than 0.2. The results for the unphysical 
range, | < ^ < 1 are included in the figures to exhibit the curious looping effect. 
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B. Discussions 



In order to study the origin of the looping effect we have systematically switched off the 
density dependences of the various coupling constants and meson masses. In so doing, we 
find that the density dependence of g'^jsi^ is the main source of the looping behavior of the 
saturation values oi E/A and p. Qualitatively such a behaviour can be understood with 
a simple model where the essential ingredient of the density dependence in the coupling 
constant ^^jvat is built in as a small modification of the saturation curve. 

Assuming that all meson masses and all coupling constants, except g'^ATAr ^-^e independent 
of pN-i the density dependence of E/A can be modelled in the following manner: 
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E/A = {E/A)o + -My - if - r^^Zili^ \{9.NN{p)/9aNN{0)f - 1 , (51) 
where, for convenience, we have introduced the quantity 

y = p/po, (52) 

The first two terms describe the density dependence oi E/A near the saturation point 
when none of the coupling constants and meson masses is density dependent. The quan- 
tity Kq is the nuclear incompressibility and {E/A)q is the minimum value of E/A in this 
situation. The last term contains the effect of the density dependence of g*^^^) calculated 



perturbatively, using, in addition, the fact that (pi — P2Y -C i.e., the average squared 
relative momentum of two Fermi sea nucleons is much smaller than the square of the cr 
meson mass. 
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FIG. 5. Plots of 2A(£;/A)/(280Mel/) vs. the saturation density p/pQ for Kq = 280MeV and 
(y^, ri) = (0.17, 0.1), (0.35, 0.2), (0.7, 0.2), and (1.0, 0.2). 

An inspection of the two graph in the left column of Fig. ^ shows that the quadratic 
form: 

gaNN{p)/gaNN{Q) = 1 + V^Vi'^ " i^) , (53) 

2|/1 

can describe these graphs quite well. The ratio gaNN^p) / QuNNiS^) peaks aX = yi and the 
peak value is 1 + ?7?/i/2. The values of the various parameters used in the illustration are 
K = 280MeV, g^NN = 8.0, = 630MeV. We consider four sets of yi and rj, namely, 
(yi, r]) = (0.17, 0.1), (0.35, 0.2), (0.7, 0.2), and (1.0, 0.2). The first two describe the ratios 
gaNN^p)/ gaNN^O) foT the QUARTIC and the PUREMASS cases, respectively, in Fig. ||. The 
last two are used to illustrate the origin of looping. 
It is convenient to recast Eq. ( pi]) into the form: 

3po gaNNjO f 
' AKq ml 



2[E/A - {E/A)o]/Ko = {y - if - 1/7^ ' {g.NN{p)/ g.NNm' - 1 • (54) 



Because of the presence of the third term in Eq. (^) the saturation density will change. 
For a given value of ^ this can be readily determined by minimizing the expression for E/A 
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with respect to y = p/po- The resulting plots of the change in the binding energy at the 
new saturation density vs. p/po are shown in Fig. ^ for the four sets of where ^ has been 
varied in the region 0. < < 0.5. 

Starting from zero density, as gaNN^p) increases with p the quantity E/A becomes more 
negative. The saturation properties reflect this through a combination of increasing satura- 
tion density, p, and decreasing E/A. Eventually, decrease of gaNN^p) beyond its peak with 
increasing p and the effect of incompressibilty take over and the saturation density stops 
increasing. If this occurs in the range < ^ < 0.5, we see looping. If the peak of gaNN^p) 
occurs at a higher density, i.e., for a higher yi, the saturation density keeps on increasing 
and we see no looping for ^ < 0.5. In the present example the critical value of yi for t] = 0.2 
appears to be slightly higher than 0.7. Thus the occurence of a peak of gaNN^p) at not too 
high a value of the effective density is needed for looping to occur for ^ < 0.5. The value 
of yi is lower for the QUARTIC case than for the PUREMASS case. This explains why the 
loop turns for the former case at a density lower than that for the latter. 

C. Concluding Remarks 

Due to the internal structure of the nucleon, we should, in general, expect that the ef- 
fective meson nucleon parameters may change in nuclear medium. Using a chiral confining 
model we have studied the density dependent changes of the meson-nucleon coupling pa- 
rameters. We have also used a simple ansatz for the density dependence of a and uj masses. 
Using the framework of a Dirac-Brueckner analysis we have found that their effect on the 
saturation properties of nuclear matter can be significant. Due to the density dependence 
of the g^NN as predicted by the chiral confining model we have found, in particular, a loop- 
ing behaviour in the nuclear matter binding energy at saturation density. The looping has 
essentially also been verified in a simple model, where it is mainly caused by the presence 
of a peak in the density dependence of the medium modified aN coupling constant at a low 
density. We should stress that the effect of the density dependence of the other quantities are 
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not negligible. However, a qualitative understanding of the looping effect can be obtained 
by paying attention to the density dependence of q^jnn alone. 

It appears that the density dependence of the coupling constants and the meson masses 
produce effects which are small but interesting. In particular the small effect tends to im- 
prove the results for nuclear matter. Prom the present study we see that the relationship 
between binding energy and saturation density may not be as universal as found in nonrel- 
ativistic studies and that more model dependence is exhibited once medium modifications 
of the basic nuclear interactions are considered. We hope that these preliminary results will 
encourage more detailed investigation of this particular variety of density dependence of the 
NN interaction. 
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